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Abstract 

In this paper we have produced different kinds of bimagic squares based on bimagic 
squares of order 8x8, 16x16, 25x25, 49x^9, etc. A different technique is applied 
to produce bimagic square of order 16x16, 25x25, 49x49, etc. The bimagic square 
of order 8x8 used is the already known in the literature. The work is neither based 
on any programming language nor on mathematical results. Just simple combina- 
tions are used to produce these bimagic squares. Moreover, in each case we have 
used consecutive numbers starting from 1. 

1 Introduction 

In this work we produced different orders of bimagic squares containing bimagic or semi- 
bimagic squares. Most of the work is based on the bimagic squares of orders 8x8, 16x16 
and 25x25. We have also produced new bimagic squares of order 49x49, 121x121. The 
bimagic square of order 8x8 is obtained long back by Pfeffermann 1891 [5]- The bimagic 
squares of order 16x16 and 25x25 are very much similar to one given in [2] as examples. 
The difference is that the bimagic square of order 16x16 appearing in [2] don't have the 
property that each sub-block of order 4x4 as a magic square, while in our case it happens. 
The bimagic square of order 25x25 is very much similar to given in [2j but we have applied 
a little different approach to produce it. The bimagic square of order 49x49 is long back 
produced by G. Tarry 1895 [Ij. Here also we applied a little different approach. We have 
produced the bimagic square of order 121x121 without knowing that is it done before. 
Based on the approach adopted in this work, we can always produce bimagic squares 
using squares of prime number such as 13^ x 13^, 17^ x 17^, etc. 

No programming language is used, just simple combinations are sufficient to produce 
whole the work. In each case, we have used consecutive numbers starting from 1. Some 
of these files are available at the authors' web-site given above. 

During construction we observe that in case of orders /c^ x A;^, k = 4, 5 and 7, in 
the previous subgroup x k^, k = 4, 5 and 7 the bimagic sums has the value in each 
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case. For more details see the table given at the end as an open problem to prove it 
mathematically. 

Before we proceed, here below are some basic definitions: 

(i) A magic square is a collection of numbers put as a square matrix, where the sum 
of element of each row, sum of element of each column and sum of each element of 
two principal diagonals have the same sum. For simplicity, let us write it as SI. 

(ii) Bimagic square is a magic square where the sum of square of each element of 
rows, columns and two principal diagonals are the same. For simplicity, let us write 
it as S2. 

(iii) Upside down, i.e., if we rotate it to 180" degree it remains the same. 

(iv) Mirror looking, i.e., if we put it in front of mirror or see from the other side of 
the glass, or see on the other side of the paper, it always remains the magic square. 

(v) Universal magic squares, i.e., magic squares having the property of upside down 
and mirror looking are considered universal magic squares. 

A good collection of multimagic squares can be seen in [1]. New upside down and 
universal magic squares can be seen in Taneja 

2 Details 

Whole the work we have divided in small parts. In the end we have given magic and 
bimagic sums in each case. 

2.1 First Part 

In this part we have presented bimagic squares of the following orders: 

16x16, 32x32, 56x56, 64x64, 72x72, 88x88, 96x96, 
104x104, 112x112, 128x128, 144x144, 176x176, 208x208, 
224x224, 256x256, 512x512, 1024x1024, 2048x2048 and 4096x4096. 

Let us divide the above bimagic squares in three small groups. 
2.1.1 First Small Group 

In this subsection we have given the following bimagic squares 

16x16, 64x64, 256x256, 1024x1024 and 4096x4096. 
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In this group we have the special property that each block of orders 16x16, 64x64, 
256x256 and 1024x1024 are also bimagic squares. Also, each block of order 4x4 is a 
magic square. In case of magic squares of order 256x256 each block of order 64x64 has 
the same bimagic sum S2. In case of magic squares of order 1024x1024 and 4096x4096 
each block of order 256x256 produces a same bimagic sum S2 for 64x64. Using the same 
procedure we can also calculate the bimagic square of order 128 x 128. Its values are given 
in the last section. 

2.1.2 Second Small Group 

In this subsection we have given the following bimagic squares 

32x32, 64x64, 96x96, 128x128, 512x512 and 2048x2048. 

This group has the special property that each block of orders 8x8, 64x64, 96x96, 
128x128, 512x512 and 2048x2048 is either a bimagic or semi-bimagic square but the 
final group is always a bimagic square. Each block of order 8x8 is always a magic square 
having the same magic sum SI in whole the order. The bimagic square of order 160x160 
can also be calculated with the same procedure, but we have calculated it in the next part 
as multiple of 16. If we calculate 128x128 according to first small group i.e., using bimagic 
square of 16x16, then the next orders 512x512 and 2048x2048 can also be calculated as 
combinations of bimagic squares of order 16x16. In this case all the bimagic squares of 
this group goes to first small group except 32x32 and 96x96. 

2.1.3 Third Small Group 

In this subsection we have given the following bimagic squares 

56x56. 72x72, 88x88, 104x104, 112x112, 
144x144, 176x176, 208x208 and 224x224. 

In case of 56x56, 72x72, 88x88 and 104x104, each block of order 8x8 is cither a 
bimagic or semi-bimagic square but the final order is always a bimagic square. Each block 
of order 8x8 is always a magic square. While, in case of 112x112, 144x144, 176x176, 
208x208 and 224x224 each block of order 16x16 is a bimagic square with the property 
that each block of order 4x4 is a magic magic square. 

2.2 Second Part 

In this part we have presented bimagic squares of the following orders: 

40x40, 80x80, 120x120, 160x160, 200x200, 240x240, 
400x400, 600x600, 800x800, 960x960, 1000x1000, 1200x1200, 
1600x1600, 2000x2000, 2400x2400, 3000x3000, 3200x3200 and 4000x4000. 

Let us divide the above bimagic squares in two small groups 
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2.2.1 First Small Group 

In this subsection we have given the following bimagic squares 

80x80, 160x160, 240x240, 400x400, 800x800, 960x960, 
1200x1200, 1600x1600, 2000x2000, 3200x3200 and 4000x4000. 

This group we have the special property that each block of orders 16x16 is a bimagic 
square with each block of order 4x4 as a magic square. In higher cases such as 400x400, 
each block of order 80x80 is also a bimagic square, etc. 

2.2.2 Second Small Group 

40x40, 120x120, 200x200, 600x600, 
1000x1000, 2400x2400 and 3000x3000, 

This group has the special property that each block of orders 8x8 is either a bimagic 
or semi- bimagic but the final group is always a bimagic square. Each block of order 8x8 
is always a magic square. In the higher case such as 200x200, each block of order 40x40 
is a bimagic or semi-bimagic, etc. 

2.3 Third Part 

In this section we have presented bimagic squares of following orders: 

25x25, 125x125 and 625x625 

Here each block of order 5x5 is a magic square with the same sum SI. Each block 
of orders 25x25 is a bimagic square. In case of 625x625, each block of order 125x125 is 
also a bimagic square with the same sum S2. 

2.4 Forth Part 

In this section we have presented bimagic squares of the following orders: 

49x49, 343x343 and 2401x2401. 

Here each block of order 7x7 is a magic square with the same sum SI. Each block 
of order 49x49 is a bimagic square. Moreover, in case of 2401x2401 each block of order 
343x343 has the same bimagic sum S2. 

2.5 Fifth Part 

In this section we have presented bimagic squares of the following orders: 

121x121 and 1331x1331. 

Here each block of order 11x11 is a magic square with the same sum SI. Each block 
of order 121x121 is a bimagic square. The next order 14641x14641 multiple of 11 with 
1331 is not calculated because of higher values. 
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3 Semi-Bimagic Squares 



Following the same approach applied to obtain bimagic squares in the above sections, we 
can still have semi-bimagic squares of orders 24x24 and 48x48 with the property that 
in case of 24x24, each block of order 8x8 is a magic square with three blocks of order 
8x8 as bimagic and six blocks as semi-bimagic. In case of 48x48, each block of order 
4x4 is a magic square. Here we have 9 blocks of bimagic squares of order 16x 16. For the 
numerical values of these semi-bimagic squares see the last section. 

4 Bimagic Squares of Magic Squares 

Tarry-Cazalas [3J in 1934 (see [Ij) gave a bimagic square of order 9x9 . Following the 
idea of Tarry-Cazalas and the approach adopted above, we obtained here bimagic squares 
of orders: 

81x81 and 729x729 . 

These two bimagic squares have the property that each block of order 9x9 is just a 
magic square, rathar than bimagic as in the other cases studied above. Also the sum of 
all 9 member in each block of order 3x3 has the same sum as of iSlgxg. If we follow the 
idea of Pfeffermann (see [I]) of bimagic square of order 9x9 and use our approach, we are 
unable to get bimagic squares of orders 81x81 and 729x729. In whole the work, this is 
the only case, where we don't have subgroups of bimagic squares. 

We can still have bimagic squares of orders 81 x81 and 729 x 729 considering just magic 
square of order 9x9. In this situation, we get bimagic squares with the property that each 
small groups of orders 3x3 and 9x9 are semi- magic squares finally giving bimagic squares 
of orders 81x81 and 729x729 . For details see the files given in authors' site [13j . 

5 Open Problem 

The bimagic square of order 16 x 16 already known in the literature don't have the property 
that each block of order 4x4 as a magic square. This we have done in this work. We have 
produced bimagic squares of order 25x25, 49x49, 121x121, etc. in a different approach 
than the one already known in the literature. Interesting the approach adopted lead us 
the following property: 
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We observe from tlie above table tliat in tlie first sub-group in eacli case the S2 is same, 
i.e., ^ ^ k'^ ^ , k = 4, 5 and 7. Now the question is to prove it mathematically? 
Moreover, if we go to higher order, for example in case of 1024x1024, the bimagic squares 
of order 64x64 has the same sum S2 in each group of 256x256, while the bimagic sums 
S2 for 256x256 give different values. 



6 Numerical Values 

Here below are the numerical values in each case. Some files of these bimagic squares can 
be downloaded at authors' web-site [I3] : 

• Bimagic square of order 16x16 

1514x4 := \SliQxl6 = 514 

SUexw := 2056 
.5216x16 := 351576 

• Semi-Bimagic square of order 24x24 

Slo^Q._i c-i„, „, = 2308 



24x24 



.51 

'5224x24 
'5224 x 24 
'5224 x 24 



6924 

2661124 - (rows and columns) 
2654292 - (diagonal - 1) 
2714116 - (diagonal - 2) 



Bimagic square of order 25x25 

'515x5 •= |'5l25x25 = 1565 
'5125x25 •= 7825 
.5225x25 := 32 6 3 25 
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• Bimagic square of order 32x32 

518x8 := i5l32 x 32 = 4 100 

5132 x 32 := 164 00 
5232x32 := 11201200 

• Bimagic square of order 40x40 

518x8 := |5l40x40 = 6404 

5140x40 ■~ 32020 
5240x40 32165340 

• Semi-Bimagic square of order 48x48 

5l4x4 := |51i6xl6 = ^5148x48 — 4610 

5l4gx48 55320 

5248x48 := 84989960 - (rows and columns) 
5248x48 := 84990120 - (diagonal - 1) 
5248x48 := 85358600 - (diagonal - 2) 

• Bimagic square of order 49x49 

5l7x7 '■— ^5149x49 — 8407 

5149x49 := 58849 
•5249x49 := 94217249 

• Bimagic square of order 56x56 

5l8x8 •= 7'S'l56 x 56 = 125 48 

5156x56 '■= 87836 
5256x56 := 183665076 

• Bimagic square of order 64x64 

5l8x8 '■— |5l64x64 — 16338 

or 

5l4x4 := |51i6xl6 = ]^'S'l64x64 — 8194 

5164x64 := 131104 
5264x64 := 358045024 

• Bimagic square of order 72x72 

518x8 := ^5172x72 = 20740 
5172x72 := 186660 
5272x72 := 645159180 
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• Bimagic square of order 80x80 

5'l4x4 := |>S'li6xi6 = ^'S'lgoxso — 12802 
5180x80 := 256040 
5280x80 := 1092522680 

• Bimagic square of order 81x81 

'S'lgxg '■= ^>S'l8ix8i = 29529 
5181x81 := 442416 
5281x81 1162527201 

• Bimagic square of order 88x88 

5l8x8 = ji^'^ssxss = 30980 
5188x88 := 340780 
5288x88 := 1759447140 

• Bimagic square of order 96x96 

5l8x8 '■= ]^5l96x96 = 36868 

5196x96 442416 
52i20xi2o := 2718351376 

• Bimagic square of order 104x104 

5l8x8 — ^51io4xio4 — 43268 
51io4xio4 := 562484 
52io4xio4 := 4056072124 

• Bimagic square of order 112x112 

5l4x4 := |51i6xl6 = ^51ii2xll2 = 25090 

51ii2xii2 := 702520 
52ii2xii2 := 5875174760 

• Bimagic square of order 120x120 

'S'lsxs '■— ^5'li2oxi2o — 57604 

5I120X120 := 864060 

52 120x120 8295264020 

• Bimagic square of order 121x121 

51iixii := Yr51i2ixi2i = 80531 
5I121X121 := 885841 
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52i2ixi2i := 8646694001 

• Bimagic square of order 125x125 

'S'l5x5 |>S'l25x25 = ^>S'li25xl25 = 39065 

-51125x125 := 97 6 6 25 
-52125x125 := 10173502625 

• Bimagic square of order 128x128 

SlsxS '■— |'S'l32x32 = ]^'S'll28xl28 — 65540 

or 

>S'l4x4 '■= I'S'lisxie = ^'S'li28xl28 = 32770 

.51128x128 := 1048640 
.52128x128 11454294720 

• Bimagic square of order 144x144 

514x4 |5'li6xl6 = ^5*1144x144 = 41474 

51i44xi44 := 1493064 
52i44xi44 := 20640614424 

• Bimagic square of order 160x160 

<5l4x4 := \SliQxi6 = ^51160x160 — 5 1202 
5I16OXI60 := 2048080 
52i60xi60 := 34954581360 

• Bimagic square of order 176x176 

514x4 |51i6xl6 — ^51i76xl76 — 61954 

51i76xi76 2725976 
52i76xi76 := 56294130376 

• Bimagic square of order 200x200 

5l8x8 •= |5l40x40 = ^51200x200 — 160004 

51200x200 := 40 00 100 
52200x200 106670666700 

• Bimagic square of order 208x208 

5l4x4 '■— |51i6xl6 — ^^51208x208 — 86530 

5 1208 x 208 4499560 
52208x208 := 129780809080 
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Bimagic square of order 224x224 

<S'l4x4 '■— I'S'liexie = ■^'S'l224x224 = 100354 

51224x224 := 56 1 98 24 
52224x224 := 187988732624 

Bimagic square of order 240x240 

1514x4 := \SliQxl6 = ^>S'l240 x 240 = 1152 02 

-51240x240 := 6912120 
-52240x240 := 265427712040 

Bimagic square of order 256x256 

<5l4x4 := |'51i6xl6 = ]^-5l64x64 — ^-51256x256 — 131074 

5164x64 := 2097184 
5264x64 := 91628066144 
-51256x560 '■= 8388736 
-52256x256 := 366512264576 

Bimagic square of order 343x343 

-517x7 '■— |-5l49x49 — ^-51343x343 — 41775 

-51343 x 343 20176975 
-52343x343 := 1582540680175 

Bimagic square of order 400x400 

>5l4x4 := |>51i6xl6 = ^-51 80 x 80 = 1^-51400x400 — 320002 

-51400x400 := 32000200 
-52400x400 := 3413365333400 

Bimagic square of order 512x512 

-5l8x8 '■— |-5l32x32 = ^-51512x512 = 1048580 

or 

<5l4x4 '■— |-51i6xl6 = ^-51i28xl28 — 1^-51512x512 — 524290 

-51512x512 := 67109120 
-52512x512 := 11728191138560 

Bimagic square of order 600x600 

-5l8x8 '■= 3^-51120x120 = ^-51 600 x 600 = 1440004 

-51600x600 := 1018000300 
-52600x600 := 25920108000100 
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• Bimagic square of order 625x625 

'S'lsxS '■— |'S'l25x25 = ^5'li25xl25 — 1^5 ^ 1625x625 — 976565 

51i25xi25 := 24414125 
52i25xi25 := 6357853190125 
51625x625 1220706625 
52625x625 := 31789265950625 

• Bimagic square of order 729x729 

5l9x9 '■— ^51729x729 — 2391489 

51729x729 := 193710609 
52729x729 := 68630571075249 

• Bimagic square of order 800x800 

5l4x4 := ^51i60xl60 = ■^51800x800 = 5120008 

51800x800 := 256000400 
52800x800 — 109226922666800 

• Bimagic square of order 960x960 

5l4x4 '■— |51i6xl6 — ^51240x240 — 2io '^"'■960x960 — 1843202 

51960 x 960 442368480 
52960x960 := 271791341568160 

• Bimagic square of order 1000x1000 

5l8x8 •= |5l40x40 = ^51 200 x 200 — 5 1 looo x 1000 — 4000004 

5I1000X1000 := 500000500 
52ioooxiooo '■— 333333833333500 

• Bimagic square of order 1024x1024 

5l4x4 '■— |51i6xl6 — ]^5l64x64 = ^51256x256 = 2^6*^^1024x1024 = 2097154 

51io24xio24 := 536871424 

52 1024x1024 := 375300505818624 

• Bimagic square of order 1200x1200 

5l4x4 '■= |51i6xl6 = ^5 1 240 x 240 — 350*^-^1200x1200 — 2880002 

5I1200X1200 := 864000600 
52i2ooxi2oo := 829440864000200 

• Bimagic square of order 1331x1331 
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'S'lllxll '■— Yr'51i2lxl21 — '511331x1331 — 9 7 4 3 591 

.511331x1331 := 1178974511 
.521331x1331 := 1392417235445951 

Bimagic square of order 1600x1600 

f5l4x4 '■— |5'li6xl6 = ^'51 80 x 80 — 1^5*1400x400 — 4^5*11600x1600 — 5120002 

5I1600XI600 := 2048000800 
52i6ooxi60o := 3495255381333600 

Bimagic square of order 2000x2000 

5*14x4 ■= 15*116x16 = ^5180x80 = i^5l4oox400 = 5^512000x2000 = 8000002 

512000x2000 := 4000001000 
522000x2000 := 106667066667000 

Bimagic square of order 2048x2048 

5l8x8 '■— |5l32x32 = ^5l5i2x512 = 216 '^■'■2048x2048 — 16777220 

or 

5l4x4 := |51i6xl6 = ^51i28xl28 — 1^851512x512 — 5^^512048 x 2048 = 8388610 

512048x2048 := 4294968320 
522048x2048 — 12009603301288960 

Bimagic square of order 2400x2400 

5l8x8 '■— ^51i20xl20 — ^51 600 x 600 = 3^0 512400x2400 — 23040004 

512400x2400 := 6912001200 

522400 x 2400 26542086912000400 

Bimagic square of order 2401x2401 

5l7x7 '■= 75I49X49 = ^51343x343 — 3^512401x2401 — 20176807 

51343x343 := 141237649 
52343x343 := 3799632717121140 
512401x2401 6920644801 
522401x2401 := 26597429019848001 

Bimagic square of order 3000x3000 

5l8x8 '■= ^51x20x120 = ^51 600 x 600 — 3^5 513000x3000 — 36000004 

5 13000x3000 := 13 5 00 1 500 
523000x3000 := 81000013500000500 

Bimagic square of order 3200x3200 
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1514x4 '■— lo'S'lieoxieo — gg'^-'-sooxsoo — 1280002 
-513200x3200 := 256000400 
-523200x3200 := 109226922666800 

• Bimagic square of order 4000x4000 

<5l4x4 '■— ]^'51i60xl60 = ■^>5l800x800 — '^''■4000x4000 — 128000008 

-514000x4000 := 32 2 00 
.524000x4000 := 341333365333334000 

• Bimagic square of order 4096x4096 

5*14x4 := 1-5116x16 = ]^-5l64x64 = ^-51256x256 = 2^6 '^■'■1024x1024 = -5 14096x4096 = 

33554434 

-514096x4096 := 34359740416 
-524096x4096 := 384307202562021376 
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